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, Abstract. In a standard method of determining electron density for soft X-ray (SXR) 

04 ' flare kernels it is necessary to assume what is the extension of a kernel along the line of 

' sight. This is a source of significant uncertainty of the obtained densities. 

Oh. In our previous paper (Bak-St§slicka and Jakimiec, 2005) we have worked out another 

method of deriving electron density, in which it is not necessary to assume what is the 
extension of a kernel along the line of sight. The point is that many flares, during their 
decay phase, evolve along the sequence of steady-state models [quasi-steady-state (QSS) 
evolution] and then the scaling law, derived for steady-state models, can be used to 
, determine the electron density. 

I The aim of the present paper is: (1) to improve the two methods of density determi- 

• ' nation, (2) to compare the densities obtained with the two methods. We have selected a 

, number of flares which showed QSS evolution during the decay phase. For these flares the 

electron density, A'', has been derived by means of standard method and with our QSS 
^ ' method. Comparison of the N values obtained with the two different methods allowed 

, us: (1) to test the obtained densities, (2) to evaluate the volume filling factor of the SXR 

C/5 ' emitting plasma. 

I Generally, we have found good agreement (no large systematic difference) between the 

values of electron density obtained with the two methods, but for some cases the values 
can differ by a factor up to 2. For most fiare kernels estimated filling factor turned out to 
^ , be about 1, near the flare maximum. 
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1. Introduction 

Loop-top flare kernels are the main sources of Soft X-ray (SXR) emission 
near flares maximum and during their decay. In order to investigate energy 
balance of the kernels it is necessary to have reliable estimates of plasma 
density inside the kernels. A "standard' method of determining the electron 



C^_' number density, consists in a few steps: 

(1) Usually, the temperature, T, is estimated from the ratio of SXR fluxes 
measured with two different filters: 

h{T) ft ^ ' 

where Fi and F2 are measured X~ray fluxes, fi{T) are response functions 
for measurements with two filters (i=l,2) and T is the estimated value 
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of the temperature. In the case of observations with the Reuven Ramaty 
High Energy Solar Spectroscopic Imager (RHESSI, Lin et al, 2002) the 
temperature is estimated from the slope of SXR spectrum. 

Flare kernels are multithermal (non- isothermal) , i.e. 

Fi = I ^{T)fi{T)dT (2) 

where '^{T) is differential emission measure (DEM). Therefore obtained 
values of T depend on instrumental response functions fi{T) and we ob- 
tain different values of T from measurements with different instruments 
(see Table I). 

(2) Next the emission measure is calculated from the formula: 

EM = Fi/fi{f) (3) 

In many cases the response function, fi{T), steeply depends on tempera- 
ture. Then random errors of determination of mean temperature, T, may 
cause significant errors of EM. To improve this standard method of EM 
determination, we recommend that observations having flat response 
functions (low \dfi/dT\ in the investigated range of temperatures) should 
be used to determine EM from Equation (3) (see Section 2). Then we 
have: 

Fi = I ip{T)fi{T)dT ip{T)dT = U EM, (4) 

i.e. we will obtain reliable EM from Equation (3), despite that the 
emitting plasma is multithermal. 

(3) The electron number density is estimated as: 



N = ^/EMJv, (5) 

where V is the volume of emitting plasma. In order to estimate the 
volume V it is necessary to assume what is the extension of emitting 
plasma along the line of sight. This introduces significant uncertainty 
into estimates of electron density from Equation (5). The uncertainty is 
high when the shape of a flare kernel is elongated (elliptical) or irregular, 
since then it is most difficult to guess what is the extension of emitting 
plasma along the line of sight. 

Additional drawback of standard method of density determination is due 
to the fact that it gives values averaged over whole volume of a kernel. In 
reality, the emitting plasma can have higher density if it occupies only 
a part of kernel volume, i.e. if its filling factor is smaller than 1. 
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In our previous paper (B§k-St§slicka and Jakimiec (2005), hereafter Paper 
I) we have proposed a new method of electron density determination in 
which (1) it is not necessary to assume what is the extension of the emitting 
plasma along the line of sight and (2) which gives actual values of the electron 
density of emitting plasma. 

Rosner, Tucker and Vaiana (1978) investigated models of steady-state coro- 
nal loops and they have found a "scaling law" of these models which can be 
written in the form: 

Nss = 2.1 X lO^rV^ (6) 

where T is the temperature at the top of a loop, L is the semilength of the 
loop and Nss is the electron number density within the loop. 
Those authors have found that Equation (6) is also valid for hot "post- 
flare loops" having temperature up to 10 MK (see their Figure 9b). This 
last finding has been explained by Jakimiec et al. (1992) who carried out 
hydrodynamic modeling of the evolution of hot flaring loops. They have 
found that if the heating rate, Enit), of a flaring loop decreases slowly during 
decay phase, the loop evolves along a sequence of steady-state models [quasi- 
steady-state (QSS) evolution] and then Equation (6) is valid during this 
evolution. Therefore this Equation (6) can be used to determine densities N 
during the QSS evolution. The aim of the present paper is: (1) to improve 
the two methods of density determination, (2) to compare the densities 
obtained with the two methods. We have selected a number of flares which 
showed QSS evolution during the decay phase. For these flares the electron 
density, A^, has been derived by means of standard method and with our 
QSS method. Comparison of the N values obtained with the two different 
methods allowed us: (1) to test the obtained densities, (2) to evaluate the 
volume filling factor of the SXR emitting plasma. 

Observations and their analysis are described in Section 2. Section 3 
contains summary and conclusions. 

2. Observations and their analysis 

Long time ago it has been shown that it is useful to display flare evolution 
in temperature-density diagnostic diagrams (Jakimiec et al. (1987, 1992)). 
In such diagrams the QSS branch of flare evolution is a straight line of 
inclination ^ = dlogT / dlogN ~ 0.5. Soft X-ray images show that size of 
loop-top sources does not change significantly during decay phase of many 
fiares {i.e their volume V ~ const). Therefore in Paper I we have found 
that simplified diagnostic diagrams (where logN is replaced by log^/EM) 
derived from GOES (Geostationary Operational Environmental Satellites) 
observations, are very useful for detection of the QSS evolution (see Fig- 
ure 1). Correct inclination, <^ ~ 0.5, of the QSS branch in these diagrams 
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indicates that temperature, TqqeSi derived from GOES observations, ad- 
equately represents mean temperature of emitting plasma during the QSS 
flare evolution and therefore we should use this temperature in Equation (6). 

In Paper I we have also found that temperatures TsxT are usually lower 
than Tqoes (SXT is the Yohkoh/Soit X-ray Telescope - see Tsuneta et at, 
1991). This is caused by the fact that SXT has narrower spectral range of 
sensitivity (1-5 keV for Be + A112 filters; 1.6-25 keV for GOES). Low values 
of TsxT result in low inclination [Q < 0.5) of the QSS branch derived with 
this temperature (see Kolomaiiski et al. (2002) and our Paper I). 

2.1. Comparison of the electron densities obtained with two 
independent methods 

Using Yohkoh/SX.T and GOES observations we have selected 20 flares which 
are limb or near-the-limb flares and their GOES diagnostic diagrams show 
clear QSS branch. For the beginning of the QSS branch (near point E in 
Figure 1) we have determined temperature, TsxT and TgoeSi and densities, 
NsxT and Nss, using the SXT images and GOES data. The flaring loop 
semilength, L, has been calculated as L = ^h, where h is the loop height 
measured in the SXR image. Obtained values of the temperatures and den- 
sities are given in Table I. Comparison of NsxT and Nss is shown in Figure 
2. We see good agreement between NsxT and Nss values. The regression 
line is: NsxT = 0.803Nss + 3.13 [this line is NsxT = 0.784iV55 + 2.21 if 
we neglect the point (30,54) which shows large deviation - see discussion 
below]. The correlation coefficient is r = 0.844. 
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Figure 1. Examples of diagnostic diagrams obtained from the GOES data. Flare evolution 
proceeds from A to F. EE is the QSS branch. Straight lines have inclination ^ = 0.5. 
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Figure 2. Comparison of Nsxt and Nss for 20 selected flares. Estimates of random errors 
are given in Table 1. 
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Table I. List of selected flares 









GOES 


GOES 






Beginning of the QSS evolution 








Date 




class 


max. 

(UT) 


Start 
(UT) 


TsXT 

(MK) 


Tgoes 
(MK) 


Nsxt 




Nss 


ANss 


k 


1. 


1991 Nov. 


04 


M3.4 


23:23 


23:47 


8.0 


11.5 


7.1 


1.4 


12.6 


2.5 


0.6 


2. 


1991 Nov. 


19 


C8.5 


09:32 


09:36 


8.0 


11.3 


9.6 


1.9 


10.4 


2.1 


0.9 


3. 


1991 Dec. 


03 


X2.2 


16:39 


16:40 


10.4 


19.1 


75.6 


15.1 


84.7 


16.9 


0.9 


4. 


1992 Jan. 


13 


M2.0 


17:34 


17:41 


8.2 


12.8 


11.9 


2.4 


15.5 


3.1 


0.8 


5. 


1992 Sep. 


09 


M3.1 


02:13 


02:15 


9.2 


14.0 


19.6 


3.9 


27.3 


5.5 


0.7 


6. 


1992 Oct. 


04 


M2.4 


22:27 


22:31 


8.1 


12.5 


14.9 


3.0 


20.3 


4.1 


0.7 


7. 


1992 Nov. 


22 


M1.6 


23:10 


23:11 


9.0 


13.3 


20.1 


4.0 


39.2 


7.8 


0.5 


8. 


1993 Jan. 


02 


C8.4 


23:52 


23:59 


7.7 


9.1 


7.2 


1.4 


9.0 


1.8 


0.8 


9. 


1993 Mar. 


23 


M2.3 


01:53 


02:11 


8.3 


10.7 


4.4 


0.9 


4.2 


0.8 


1.0 


10. 


1994 Jan. 


16 


M6.1 


23:25 


23:35 


9.0 


12.3 


30.4 


6.1 


20.4 


4.1 


1.5 


11. 


1998 Aug. 


18 


X2.8 


08:24 


08:31 


10.9 


15.6 


40.5 


8.1 


30.1 


6.0 


1.3 


12. 


1998 Aug. 


18 


X4.5 


22:15 


22:27 


11.1 


17.1 


54.1 


10.8 


30.0 


6.0 


1.8 


13. 


1998 Aug. 


19 


M3.0 


14:26 


14:50 


8.0 


10.3 


14.9 


3.0 


6.8 


1.4 


2.2 


14. 


1998 Nov. 


22 


X3.7 


06:42 


06:45 


11.2 


16.5 


38.1 


7.6 


52.7 


10.5 


0.7 


15. 


1998 Nov. 


22 


X2.5 


16:23 


16:28 


12.3 


16.3 


35.5 


7.1 


34.5 


6.9 


1.0 


16. 


1998 Nov. 


23 


X2.2 


06:44 


06:57 


11.3 


15.3 


27.4 


5.5 


26.5 


7.3 


1.0 


17. 


1999 Aug. 


04 


M6.0 


05:57 


06:09 


8.5 


12.9 


17.4 


3.5 


22.0 


4.4 


0.8 


18. 


2000 Jan. 


12 


Ml.l 


20:49 


20:57 


8.1 


10.6 


8.1 


1.6 


11.4 


2.3 


0.7 


19. 


2000 Jun. 


23 


M3.0 


14:31 


14:42 


9.0 


12.3 


14.0 


2.8 


10.1 


2.0 


1.4 


20. 


2001 Oct. 


29 


M1.3 


01:59 


02:01 


8.3 


13.2 


18.5 


3.7 


38.4 


7.7 


0.5 


Values of Nsxt and Nss are in IO^^ctt; 




Nsxt /Nss, AN 


are 20% random errors (see text) 






10 



20 



T [MK] 



30 



Figure 3. EM for various values of temperature. The calculations have been carried out 
for a constant value of the SXR flux, Ft, measured with the A112 filter {Ft — 1.5 x 10^ 
DN). 



The good agreement between NsxT and Nss values shows that both 
methods of density determination give rehable values of the density (no 
large systematic errors in their calibration). 

It is surprising that we obtain agreement between the densities despite of sig- 
nificant systematic differences between the temperatures TgxT and Tgoes- 
Explanation of this paradox is the following: The response function, fi(T) 
for the Yohkoh/SXT with A112 filter changes only a little in wide range of 
temperature, T w 8—40 MK [see curve / in Figure 9 in Tsuneta et al. (1991)]. 
Therefore, the emission measure, calculated according to Equation (3), only 
weakly depends on temperature T. To illustrate this effect, we have assumed 
a constant value of the X-ray flux, Fj, and calculated EM for various values 
of temperature (Figure 3). For A112 filter the dependence of EM on T is 
weak in the range of flare temperatures. Dependence of the electron density, 
^SXT, on temperature is even weaker, since NsxT ~ V EM. This explains 
why we obtain correct values of NsxT, despite of the fact that TgxT and 
\dTsxT /dt\ are systematically too low. 

Figure 3 displays function 1.5x 1Q^[DN]/ fi{T) for SXT observations with 
A112 filter. We have calculated values of this function for TsxT values which 
are presented in Table 1. From these values of the function we obtained 
its mean value and random (r.m.s.) error. This gives the following form of 
Equation (3): 
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EM[l0^^cm-^] = (6.48 ± 0.27) x 10^^ Fi[DN] (7) 

which does not depend on estimated value of mean temperature T, provided 
that T > 8 MK (DN is Digital Number in Yohkoh/SXT observations). 
However, the digital coefficient in Equation (7) can slowly change with time 
due to slow change in calibration of the instrument. 

In Table I we see that most flares in our sample were strong ones (GOES 
class M and X). Electron number density, measured near SXR flare maxi- 
mum, for most of the flares falls within the limits of (4 — 50) xlO^*^ cm 
But for one flare [No. 3 in Table I, point (85,76) in Figure 2] the density was 
extremely high, N ^ 8 xlO^^ cm This high density has been obtained 
with both methods of density determination. SXR image of this flare is 
shown in Figure 4. We see that it was a low flaring loop with very bright, 
compact loop-top kernel (at the beginning of the QSS evolution the diameter 
of the kernel was about 6000 km). This flare was also very strong in HXRs 
during its impulsive phase (about 1300 counts/s/SC in the Yohkoh/HXT 
23-33 keV light curves, SC is a subcollimator). 

Main source of random errors (deviations from the regression line in Fig- 
ure 2) in the standard method of density determination from SXR images, 
is due to the fact that we should assume what is extension of a SXR kernel 
along the line of sight. Main source of random errors in the QSS method 
of density determination is due to uncertainty in estimates of the flaring- 
loop semilength, L. Estimates of L are usually based on assumption that 
the shape of the loop is semicircular and (in the case of limb flares) that 
the loop lies in the plane of image. Deviations from these assumptions are 
the main source of random errors in the QSS method of electron density 
determination. 

Next we have calculated values of the ratio k = Nsxt/^ss- The values of 
the ratio k are given in Table I and their distribution is shown in Figure 5. We 
see that for two flares (No. 12 and 13) the values of k are highest {k = 1.8 and 
2.2). SXR images of these two flares are shown in Figure 6. We see that their 
SXR kernels are elongated. The volumes, V, of such kernels were calculated 
under assumption that they have "sausage-like" shape i.e. that they are 
ellipsoids with semiaxes axbxb. The high values of the ratio k suggest that 
actual volume of the SXR emitting plasma is significantly larger, i. e. the 
third semiaxis is larger than b. This indicates that estimates of NsxT for 
such elongated kernels can be loaded with highest random errors. Therefore, 
we have excluded these two large values of k from calculation of the mean 
value k. We have obtained: 

k = 0.88, s = ±0.28 (8) 

where s is random (r. m. s.) error of single estimate of k. 
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16:00 16:20 16:40 17:00 17:20 17:40 24.4 24.5 24.6 24.7 24.8 24.9 25.0 
Stort Time (03-Dec-91 15:00:00) logVEM 

Yohkoh/SXT Bel 19 3-Dec-91 16:40 UT 




-940 -920 -900 -880 -860 

X (arcsecs) 

Figure 4- GOES light curve, diagnostic diagram and sample SXT(Bell9) image (inverse 
intensity scale in colour) of December 3, 1991 flare (solid line shows solar limb). 



In Figure 5 we see that fc-values are concentrated around this mean value. 
The value of k which is close to one, confirms the good agreement between 
^SXT and Nss values. We see that random [r.m.s.) error, s, of A;-values is 
about 30%. Assuming that contribution of random errors of NgxT and Nss 
to this random error of k are similar in size, we conclude that random error 
of a single determination of NsxT or Nss is about 30%/\/2 « 20%. 

The values of NsxT, obtained with the standard method, are mean 
electron number densities averaged over the SXR-kernel volume. But the 
QSS method of density determination provides actual density of the emit- 
ting plasma. Therefore, the ratio, k = Nsxt/Nss, provides an estimate 
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Figure 5. Histogram of the ratio k — Nsxt /Nss for 20 analysed flares 

of the fraction of kernel's volume which is filled with the SXR emitting 
plasma ("filling factor"). In Figure 5 we see that A;- values are concentrated 
around k ~ 0.9 which means that, typically, near the flare maximum the 
kernel's volume is nearly entirely filled with the hot SXR emitting plasma. 



Yohkoh/SXT Bel 19 18-Aug-98 08:27 UT Yohkoh/SXT Bel 19 19-Aug-98 14:52 UT 




-840 -820 -800 -780 -760 -840 -820 -800 -780 -760 

X (arcsecs) X (orcsecs) 



Figure 6. Sample SXT(Bell9) images (inverse intensity scale in colour) of August 18, 
1998 {left) and August 19, 1998 {right) flares 
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Remaining, (1 — k), fraction of the volume can be filled with plasma of 
lower temperature, T < 4 MK, whose contribution to the investigated SXR 
emission is marginal. Values k > 1 are due to random errors in both methods 
of density determination. 



3. Summary and conclusions 

We have found that simplified diagnostic diagrams (T vs. VeM), derived 
from GOES observations, are adequate to identify quasi-steady-state (QSS) 
phase of flare evolution. Correct inclination, ~ 0.5, of the QSS branches in 
the diagnostic diagrams shows that the temperature, Tgoes, derived from 
the GOES observations adequately represents mean temperature of the hot, 
SXR emitting plasma. 

For the QSS branch of flare evolution we were able to determine electron 
number density, N, for SXR loop-top flare kernels, with two independent 
methods which are described in Sections 1 and 2. Good agreement of the 
A^- values obtained with the two methods confirms reliability of the densities 
derived from Yohkoh/SX.T images. In Section 2 we have explained why we 
obtain correct NsxT electron densities, despite of the fact that Yohkoh/SXT 
temperatures are systematically too low. 

We have obtained that electron density near the SXR flare maximum falls 
within the limits of (4 — 50)xl0^'^ cm~^, but in one case of strong, compact 
flare it was extremely high, ~ 8 xlO^^ cm~^. 

Densities, NsxT, derived from SXR images, are mean electron densities 
averaged over SXR flare kernels. But the densities, Nss, obtained with the 
second (QSS) method, provide actual densities of emitting plasma. Therefore 
the ratio Nsxt/^SS provides an estimate of the filling factor of the kernel 
with hot, SXR emitting plasma. Our analysis of obtained k values shows 
that, typically, flare kernels near SXR flare maximum are nearly entirely 
filled with the hot, SXR emitting plasma. 
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